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Abstract. We describe a new tool for the analysis of the 
UV to the sub-millimeter (sub-mm) spectral energy dis- 
tribution (SED) of spiral galaxies. We use a consistent 
treatment of grain heating and emission, solve the radi- 
ation transfer problem for a finite disk and bulge, and 
self-consistently calculate the stochastic heating of grains 
placed in the resulting radiation field. 

We use this tool to analyse the well-studied nearby 
edge-on spiral galaxy NGC891. First we investigate 
whether the old stellar population in NGC891, along with 
a reasonable assumption about the young stellar popu- 
lation, can account for the heating of the dust and the 
observed far-infrared and sub-mm emission. The dust dis- 
tribution is taken from the model of Xilouris et al. (1999), 
who used only optical and near-infrared observations to 
determine it. We have found that such a simple model 
cannot reproduce the SED of NGC891, especially in the 
sub-mm range. It underestimates by a factor of 2-4 the 
observed sub-mm flux. 

A number of possible explanations exist for the missing 
sub-mm flux. We investigate a few of them and demon- 
strate that one can reproduce the observed SED in the 
far-infrared and the sub-mm quite well, as well as the ob- 
served radial profile at 850 /an. 

For the models calculated we give the relative propor- 
tion of the dust radiation powered by the old and young 
stellar populations as a function of FIR/sub-mm wave- 
length. In all models we find that the dust is predomi- 
nantly heated by the young stellar population. 
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1. Introduction 

Dust grains can be considered as test particles for the 
intrinsic radiation field in galaxies. Observations of their 
emission in the infrared (IR), combined with optical and 
ultraviolet (UV) observations of the light from stars, at- 
tenuated by the grains, should, in principle, strongly con- 
strain the intrinsic distribution of stellar luminosity and 
dust in galaxies. This would address the fundamental 
question of optical thickness in galactic disks and allow 
evaluation of optical and IR observational data in terms 
of physical (e.g., star formation history) rather than em- 
pirical (e.g., colours) parameters. Ultimately, such analy- 
sis techniques will be necessary to interpret the evolution 
of the UV - far- infrared (FIR)/sub-millimeter (sub-mm) 
spectral energy distribution (SED) of galaxies and the IR 
background radiation over cosmological timescales. 

However, even for galaxies in the local universe, techni- 
cal difficulties are considerable from both theoretical and 
observational stand points. From the theoretical point of 
view, it is necessary to solve the radiation transfer problem 
for a finite disk, and self consistently calculate the stochas- 
tic heating of grains placed in the resulting radiation field 
using realistic grain models. To that must be added the 
problem that disks are fundamentally inhomogeneous on 
different scales. From the observational point of view, it 
is only now becoming possible to obtain resolved images 
covering the mid-IR to sub-mm spectral range for nearby 
galaxies. 

Because of these difficulties, the origin of the FIR emis- 
sion in local universe spiral galaxies is still a subject of de- 
bate. Based on models of dust heating in various environ- 
ments (e.g., Mezger et al. 1982; Cox et al. 1986), and from 
comparison of IRAS data with optical properties of galax- 
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ies, several studies concluded that a substantial fraction 
of the FIR emission from galaxies could be due to dust 
heated by the diffuse interstellar radiation field (ISRF), 
the remaining part being due to localised sources around 
OB stars in star-forming regions (e.g., de Jong et al. 1984; 
Walterbos & Schwering 1987). Some later papers, how- 
ever, suggested that the correlation between Ha and FIR 
emission implies that OB stars by themselves can account 
for all the dust heating, so that there is no need for an 
additional contribution to heating by the ISRF (e.g., De- 
vereux & Young 1990). By contrast, there have been a 
series of papers (e.g., Xu et al. 1994) which concluded 
that the non-ionising UV accounts for most of the en- 
ergy absorbed by dust. They showed that the FIR/radio 
correlation can be separated into a warm-FIR/thermal- 
radio correlation which is due to massive ionising stars 
(> 20 M Q ), and a cool-FIR/nonthermal-radio correlation, 
principally due to intermediate mass stars (5 — 20 M Q ) - 
the supernova progenitors. A heating component from the 
old stellar population was also needed to account for the 
non-linearity of the correlation, but this accounts only for 
about a third of the heating. 

In this paper we describe a new method for the analysis 
of the radiation energy budget in edge-on spiral galaxies, 
as part of an effort to develop a theoretical tool for the 
interpretation of the SED in spiral galaxies. The method 
presented here is based on some knowledge of the primary 
radiation field (i.e., from stars) in spiral galaxies and the 
dust distribution in them. For edge-on spirals, it has been 
demonstrated (Xilouris et al. 1997; 1998; 1999) that it 
is possible to determine the distribution of older stars 
and associated dust from analysis of optical and near- 
infrared (NIR) images. Thus, the diffuse radiation field in 
the galaxy from these optical bands can be computed quite 
accurately. The diffuse radiation field in the UV is gener- 
ally not so well constrained by the UV data alone, due to 
the larger optical depth of disks in the UV. This younger 
stellar component is parameterised in terms of a recent 
star-formation rate (SFR) and it can be observationally 
constrained in the case that a direct, optically thin indi- 
cator, such as radio free-free emission, is available. Solving 
the radiation transfer problem within the disk and bulge, 
and calculating the heating of grains we can determine the 
diffuse FIR/sub-mm emission. Potentially, further contri- 
butions arising from emission from localised sources need 
to be added to the diffuse component. Then, because the 
dust is optically thin in the FIR, the predicted emission 
can be compared with observations. 

We apply the above method to the well-known edge-on 
spiral galaxy NGC891, for which we assume a distance of 
9.5 Mpc (van der Kruit & Searle 1981). This is one of the 
most extensively observed edge-on galaxies in the nearby 
universe, which makes it ideal for a verification of our 
modelling technique. The analysis of a sample of 6 more 
edge-on galaxies is in progress (Misiriotis et al. 2000a). 
Modelling edge-on spiral galaxies has several advantages, 



mainly when investigating them in the optical band. One 
advantage is that, in this view of a galaxy, one can easily 
separate the three main components of the galaxy (i.e., 
the stellar disk, the dust and the bulge). Another advan- 
tage is that the dust is very prominently seen in the dust 
lane, and thus its scalelength and scaleheight can be bet- 
ter constrained. A third advantage is that many details of 
a galaxy that are evident when the galaxy is seen face- 
on (e.g., spiral arms), are smeared out to a large degree 
when the galaxy is seen edge-on (Misiriotis et al. 2000b). 
Thus, a simple model with relatively few parameters can 
be used for the distribution of stars and dust in the galaxy. 
However, the third advantage comes with disadvantages, 
especially when trying to model the galaxy in all wave- 
length ranges, including UV and FIR/sub-mm. Thus, in 
edge-on galaxies it is very difficult to see localised sources 
(i.e., HII regions), in which the radiation can be locally 
absorbed and thus not contribute to the diffuse radiation 
field. Also, if the galaxy has a thin (young) stellar/dust 
disk, highly obscured by the dust lane in the plane of 
the galaxy, then this disk cannot be inferred from obser- 
vations in the optical/NIR spectral range. In passing we 
mention that throughout this paper we will use the terms 
"thin/thick disks" to describe their scaleheights, and not 
their optical thickness. 

The goal of this paper is to find the simplest phe- 
nomenological description of NGC891 that can adequately 
account for the observed optical to sub- mm SED. In par- 
ticular we investigate what fraction of the FIR is ac- 
counted for by the Xilouris et al. (1999) model, namely 
by a disk of dust, heated by a diffuse, smooth, optically 
emitting, old stellar disk and bulge, supplemented by a dif- 
fuse, smooth, UV emitting disk of newly-formed stars. We 
refer to this model as the "standard model" . We will show 
that this model fails to reproduce the observed SED of the 
galaxy at FIR/sub-mm wavelengths. We investigate then 
different possibilities to explain the origin of the missing 
FIR/sub-mm flux. 

Recently there have been several works modelling the 
SED of galaxies from the UV to the sub-mm (Silva et 
al. 1998, Devriendt et al. 1999). Their method consists of 
using models for photometric and/or spectrophotometric 
and chemical evolution of galaxies in order to fit the ob- 
served SED. For example Silva et al. (1998) used a chem- 
ical evolution code to follow the SFR, the gas fraction, 
and the metallicity of the galaxies. In their approach the 
parameters describing the star formation history as well 
as the geometrical parameters of the intrinsic distribution 
of stars and dust are left as fitting parameters. While this 
approach was very successful in fitting the SED of dif- 
ferent types of galaxies, it implies the use of many free 
parameters. Furthermore, a more detailed description of 
the geometry of the galaxy is indispensable for an accurate 
determination of the dust content in individual cases. 

This paper is organised as follows: In Sect. 2 we de- 
scribe our "standard model" for the diffuse dust and ra- 
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diation field, including the optical properties of the dust 
grains, and the derivation of the diffuse radiation field 
powered by the old and the young stellar populations. In 
Sect. 3 we give the results of our calculations, in terms 
of integrated SED for the "standard model" and in Sect. 
4 we supplement our model with the contribution of lo- 
calised sources. Since the "standard model" fails to repro- 
duce the observed SED, we discuss in Sect. 5 the origin 
of the FIR emission, considering four possibilities to pro- 
duce the missing IR component. From these, the two dust 
disk model, together with the clumpy scenario constitute 
the simplest solutions we have been able to identify to 
describe the origin of the missing infrared component. In 
Sect. 6 we give our results in terms of radial profiles and in 
Sect. 7 we discuss the contribution of radiation at different 
wavelengths in heating the dust. In Sect. 8 we summarise 
our results and give our conclusions. 

2. The model 

2.1. The properties of the dust grains 

We consider only the spectral range A > 40 /im as our dust 
model does not include Polycyclic Aromatic Hydrocarbon 
(PAH) molecules, which are important for A < 15 /im and 
may even contribute at A = 25 /im, a point that is still very 
uncertain. For example, Desert et al. (1990) proposed a 
dust model in which 48% of the IR emission at 25 /im (for 
dust in the solar neighbourhood) is emitted by the PAH 
molecules as continuum emission. Draine and Anderson 
(1985) proposed that the IR emission at 25 /im and shorter 
wavelengths can be explained by enhancing the number 
of small grains. Our restriction to the spectral range A > 
60 /im allows us to avoid the uncertainties introduced by 
a different dust component (PAHs or modified grain size 
distribution). 

In characterising the dust properties we thus consid- 
ered spherical "astronomical grains" for two component 
materials - graphite and silicate - of the Mathis et al. 
(1977; hereafter MRN) interstellar grain model. MRN 
have proposed that the observed interstellar extinction of 
star-light may be produced by a mixture of graphite and 
silicate particles with a simple power-law distribution of 
sizes 



dNi = Ni(a) da = Nid k da 



(1) 



where dNi is the number of grains of type i with radii in 
the interval [a, a + da], and k = 3.5. We consider the lower 
and the upper cutoffs to the size distribution to be respec- 
tively a m i n = 10 A and a max = 0.25 /im. The graphite and 
silicate abundances were taken from Draine & Lee (1984, 
hereafter DL), and are only slightly different from those 
found by MRN; DL proposed a mixture of 53% silicates 
(A^si) and 47% graphites (Acraphite), which were chosen 



to fit the extinction curve in our Galaxy, and which we 
also adopted for NGC891. 

The absorption efficiencies Q a bs(A,a), the scattering 
efficiencies Q SC at(A,a) and the scattering phase function 
g(X, a) were taken from Laor & Draine (1993). The ab- 
sorption coefficient K a b s , the scattering coefficient K scat as 
well as the averaged anisotropy parameter of the Henyey- 
Greenstein scattering phase function g can then be derived 
by integrating Q abs (A, a), Q sca t(A,a), and g(X,a), respec- 
tively, over the size distribution 

^2 A; a~ k cta7ra 2 Q abs (A, a) (2) 

-mi„ l 

Kscat(A) — / A; a~ k daira 2 Q sca t(A,a) (3) 



/' 



^A;a k daTra 2 g(X,a)Q scat (X,a) 

.'/(A) = ''"""a,! (4) 

^Nia k da-Ka 2 Q scat {X, a) 



with the extinction coefficient n cxt = n a b s + «scat • 

As the 60 — 100 /im data are expected to contain a sig- 
nificant amount of emission from small grains not in equi- 
librium with the radiation field, it is necessary to model 
the stochastic emission from the dust. For this calcula- 
tion we adopted the heat capacities of silicate grains from 
Guhathakurta & Draine (1989), which were derived as a 
fit to experimental results for Si02 and obsidian at tem- 
peratures 10 < T < 300 K (Leger et al. 1985), with a 
simple extrapolation for T > 300 K. For graphite grains 
the heat capacities were taken from Dwek (1986). 

2.2. Derivation of the diffuse optical and UV radiation 
field. 

Our first step is to estimate the diffuse radiation field at 
any point in the galaxy. The most straightforward way to 
do this is to solve the radiative transfer equation (RTE) for 
a given distribution of emitters and absorbers of radiation. 
In this section we will consider the distribution of emit- 
ters and absorbers as derived from optical wavelengths by 
Xilouris et al. (1998,1999), while at the UV wavelengths 
we will fix the parameters describing the geometry of the 
emitters and use population synthesis models to parame- 
terise the amplitude of the UV radiation. 

2.2.1. Radiation transfer 

In this subsection we will describe our solution of the RTE. 
Our method is based on that of Kylafis & Bahcall (1987). 
This method will be briefly described here for the sake of 
completeness. Let us introduce the RTE (Mihalas 1978) 



4 



C.C.Popescu et al. 



for a time independent radiation field along a line of sight, 
at a fixed wavelength 



dl(s, h) 
ds 



-K ext I(s, n) + r] (s, h) + r] s (s, h) 



(5) 



where I{s, h) is the specific intensity of the radiation field 
at a point s on the line of sight (defined by n), n ext is the 
extinction coefficient for the wavelength considered, ds is 
an infinitesimal length clement along the line of sight, r)o 
is the emissivity along the line of sight due to primary 
sources (stars) and rj s is the emissivity due to scattering 
into the line of sight. For coherent scattering, rj s is given 

by 



i) s (s,n)=WK el( J I(s,h')p(h,n) 



dVL' 
4ir 



(6) 



where p(n, h') is the Henyey-Greenstein phase func- 
tion (Henyey & Greenstein 1941), parameterised by the 
anisotropy parameter g, and lo is the scattering albedo 
for the wavelength considered, defined as lo = n S cat/n e xt- 
Let I = Iq + I\ + ... + I n + where Iq is the intensity 
of photons that come directly from the primary sources, 
I\ the intensity of photons that have been scattered once 
after leaving the pnirieiry sources £ md /„ the intensity of 
photons that have suffered n scatterings. The intensities 
/„ satisfy the following equations (Henyey 1937) 



dl n (s,h) 
ds 



= -K ext I n (s,n) +r] n (s,n) 



(7) 



where r] n (s,h) is the emissivity due to photons scattered 
n times, and is given by 



rj n (s,n) = ujK ext J I n _i(s,n')p(h,n')^- 



(8) 



Eq. (7) can be solved for as high order of n as allowed 
by the given hardware and/or software. In our method, 
we calculate Jo and 1\ and approximate the higher orders 
according to Kylafis & Bahcall (1987). Therefore, the so- 
lutions are given by 



io,i(s,n)= / ds' r/o,i(s', n) exp[-r(s, s')] (9) 
Jo 



while for n > 1, 
h 



In+l — In 



(10) 



where t(s, s') is the optical depth between s and s' . Given 
the intensity of the radiation field at a point, the energy 
density can be calculated by (Mihalas 1978) 



i(x,y,z) = c 1 j I(x,y,z,h) 



(11) 



2.2.2. Optical radiation field and dust distribution 

The determination of the intrinsic stellar content of a 
galaxy is often difficult because of the uncertainty intro- 
duced by dust obscuration. However, in the case of edge- 
on galaxies, dust obscuration can be used to our benefit 
to determine the spatial distribution of dust. In an edge- 
on galaxy, dust becomes visible as a very prominent dust 
lane along the major axis of the galaxy. The effects of the 
dust are strong enough to allow the existence of a fairly 
constrained model for the distribution of old stars and as- 
sociated dust (Xilouris et al. 1997; 1998;1999; Kuchinski 
et al. 1998; Ohta & Kodaira 1995). This does not include, 
however, any stellar population with a scaleheight small 
enough to be effectively completely hidden by the dust. 
Such a stellar population will be discussed in Sect. 2.2.3. 

NGC891 was modeled (Xilouris et al. 1998,1999) by 
fitting an artificial image produced by radiative transfer 
on observations in an attempt to conclude on the dust 
morphology. In this section we will adopt the latest model 
for NGC891 (Xilouris et al. 1999). In this model, the emis- 
sivity is described by an exponential disk and a de Vau- 
couleurs bulge, while the dust resides in a pure exponential 
disk. The parameters describing the stars and the dust are 
determined in 5 wavebands (B, V, I, J, K). 

For the stellar emissivity the following formula is used 



r] (R, z) = L s exp 



h. 



+ L fe exp(-7.67B 1/4 )B- 7/8 



(12) 



In this expression the first part describes an exponential 
disk, and the second part describes the bulge, which in 
projection is the well-known i? 1 / 4 -law (Christensen 1990). 
Here R and z are the cylindrical coordinates, L s is the stel- 
lar emissivity per unit volume per steradian at the center 
of the disk and h s and z s are the scalelength and scale- 
height respectively of the stars in the disk. For the bulge, 
Li, is a normalisation constant, while B is defined by 



B = 



^JR 2 + z 2 {a/by 

Re 



(13) 



with R e being the effective radius of the bulge and a and 
b being the semi-major and semi-minor axis respectively 
of the bulge. In this paper we will refer to the exponential 
disk described above as to the old stellar disk, in order to 
differentiate it from the young stellar disk emitting in UV 
(Sect. 2.2.3). 

For the dust distribution a similar expression as that 
adopted for the stellar distribution in the disk is used, 
namely 



K ext (X, i?, z) = K ext (\, 0, 0) exp ( - — 

h-d z d 



(14) 



where K ext {\, 0,0) is the extinction coefficient at the cen- 
ter of the disk for the wavelength considered, and hd and 
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Table 1. Parameters of NGC891 derived by Xilouris et 
al. (1999). 

Par am. 
Ls [- 
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Zd are the scalelength and scaleheight respectively of the 
dust. The central optical depth of the model galaxy seen 
face-on is 



r^(A) = 2 /w (A, 0,0) 2d 



(15) 



The parameters determined for NGC891 are presented in 
Table 1. The total amount of dust derived from this model 
is M dust = 5.6 x 10 7 M Q (Xilouris et al. (1999). 

The extinction coefficients K ext {\, 0, 0) were derived as 
fitting parameters (Xilouris et al. 1999), separately for 
each of the five wavebands used for modelling the galaxy. 
They were then used to calculate the absorption of the 
optical light by dust, via the radiation transfer procedure. 
Since in calculating the emission of the dust we used the 
absorption efficiencies (Q a bs) taken from Laor & Draine 
(1993), as presented in Sect. 2.1, we compare in Fig. 1 the 
extinction coefficients derived from the fitting routine with 
the theoretical ones derived from the extinction efficien- 
cies (Eqs. [2] and [3]). The result of such a comparison 
shows that indeed the extinction derived from observa- 
tions (square symbols) is very close to the theoretical one 
used to calculate the IR emission (solid line) . This ensures 
that the treatment of the absorption and of the emission 
is consistent. In Fig. 1 we also give the wavelength depen- 
dence of the absorption and scattering coefficients, plotted 
with dashed and dotted lines, respectively. 

2.2.3. UV radiation field 

The UV radiation field cannot be derived directly from ob- 
servations. For edge-on galaxies it is extremely difficult to 
infer something about the young stellar disk, which, due 
to it's small scaleheight, is presumably highly obscured 
by line-of sight dust. Any young stellar population whose 
scaleheight is smaller than the scaleheight of the diffuse 
dust will be totally unprobed by UV, optical and NIR data 
on an edge-on galaxy. Thus for an edge-on galaxy mid to 
far-IR emission would be the primary observational signa- 
ture for any young massive star population having a scale- 
height of 100 pc or less. This stellar population would (we 
assume luminosity is dominated by stars) produce almost 
all the non-ionising and ionising UV from the galactic disk. 



10" 



10" 




Fig. 1. The wavelength dependence of the extinction 
(solid line), absorption (dashed line) and scattering (dot- 
ted line) coefficients derived from the extinction, absorp- 
tion and scattering efficiencies taken from Laor & Draine 
(1993). The symbols represent the extinction coefficients 
derived as fitting parameters by Xilouris et al. (1999). 



We thus consider the UV radiation field as a parame- 
ter and we parameterise the UV luminosity in terms of a 
recent SFR, based on the population synthesis models of 
Mateu & Bruzual (2000). We consider Z = Z , a Salpeter 
IMF, M up = 100 Mq and r = 5 Gyr. This gives, for ex- 
ample, SFR = 8.55 x 10 - 28 Iw [erg/s/Hz] for A = 912 A, 
SFR = 3.74 x I0~ i3 L uv [erg/s] for the integrated ion- 
ising UV shortwards of 912 A, SFR = 1.18 x 10" 28 L uv 
[erg/s/Hz] at A = 1500 A, etc. The diffuse UV radiation 
derived from the calibration above was confined in a thin 
disk, with a scaleheight z = 90 pc (close to that of the 
Milky Way, Mihalas & Binney 1981) and the same scale- 
length as the blue disk. 

The UV luminosity could, in principle, also be indi- 
rectly inferred from different indicators of the SFR, like 
the free-free emission or Ha emission. There are however 
large uncertainties in deriving these quantities, especially 
in the case of NGC891. The first indicator, the radio free- 
free emission, has the advantage of being optically thin. 
Niklas et al. (1997) found that the total radio emission 
of NGC891 is dominated by synchrotron emission, while 
the free-free emission constitutes a thermal fraction of just 
0.05 at an observing frequency of 10 GHz. However, the 
radio morphology of NGC891 is rather complicated, with 
more morphological components having different spectral 
indices. Allen et al. (1978) suggested the existence of two 
components for the radio continuum emission: a highly 
flattened thin disk component coinciding with the equa- 
torial plane of the galaxy and a thick disk or halo com- 
ponent (see also Hummel et al. 1991). If this is the case, 
then the free-free emission estimated from a spectral de- 
composition of the integrated fluxes of the galaxy (as done 
by Niklas et al. 1997) could have been seriously underes- 
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Fig. 2. Examples of temperature distribution for dust grain sizes (10A, 20 A, 50 A, ...) embedded in the diffuse 
radiation field, for a solution with the UV radiation field in the diffuse component corresponding to a SFR = 3M©/yr: 
a) silicate grains in the center of the galaxy (R = Okpc, z — Okpc); b) graphite grains in the center of the galaxy; c) 
silicate grains at the edge of the galactic disk (R = f5kpc, z = Okpc); d) graphite grains at the edge of the galactic 
disk. 



timated. We conclude that the free-free emission does not 
constitute an accurate constraint of our model. 

The second indicator of the SFR is the Ha emission. It 
is known that correcting Ha data for internal extinction is 
a very uncertain step, since the extinction contains a large- 
scale component plus extinction due to localised sources, 
which are difficult to account for in edge-on disks. Hoopcs 
et al. (1999) found that in NGC891, 83 - 86% of the ob- 
served Ha emission is in the diffuse component associated 
with the halo of the galaxy. This would imply that most 
of the emission from the disk is highly obscured, making 
unreliable any derivation of the SFR in the disk. 

2. 3. Calculation of the IR spectrum from grains embedded 
in the diffuse optical and UV radiation field 

The dust heating due to photons from the diffuse inter- 
stellar radiation field was calculated following the method 



of Guhathakurta & Draine (1989). This method derives 
the temperature distribution P(a, T) of various grain radii 
a as a function of dust temperature T. Due to stochastic 
heating, small grains undergo significant fluctuations from 
the equilibrium temperature, while larger grains have nar- 
rower probability functions, eventually approaching delta 
functions. We calculated the temperature distribution at 
any point in the galaxy, based on the energy density of 
the radiation field derived in the previous subsections. The 
spectral shape of the UV-NIR energy density and its mag- 
nitude will obviously determine the shape of the tempera- 
ture distributions. The grain sizes used in the calculations 
were chosen between a m i n = 10 A and a max — 0.25 /zm, 
for a step Alog a = 0.05. In Fig. 2 we give the temperature 
distribution as dp/dlogT for some silicate (Fig. 2a,c) and 
graphite (Fig. 2b, d) grain sizes embedded in the diffuse ra- 
diation field, for a solution with the UV radiation field in 
the diffuse component corresponding to a SFR = 3 M Q / yr. 
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Fig. 3. The infrared spectrum from dust embedded in the 
diffuse radiation field, for a solution with the UV radiation 
field in the diffuse component corresponding to a SFR 
= 3M0/yr. The figure illustrates the infrared spectrum 
obtained for two extreme locations in the galaxy: in the 
center of the galaxy (R = Okpc, z — Okpc) - plotted with 
solid line and at the edge of the galactic disk (R = 15 kpc, 
z = kpc) - plotted with dashed-dotted line. 

Fig. 2 also illustrates the dependence of the grain temper- 
ature on the position in the disk. We give two examples, 
for two extreme positions in the plane of the disk: in the 
center of the disk (Fig. 2a, b) and at the edge of the disk 
(Fig. 2c, d). At the edge of the disk the probability distri- 
bution converges towards lower equilibrium temperatures, 
as expected due to the decrease in the intensity of the ra- 
diation field. We also note the relatively large difference 
in the equilibrium temperature between ~ 30 K for the 
central region and ~ 15 K for the outer disk. This is inter- 
esting in the context of the results of Alton et al. (1998) 
who fitted the SED of NGC891 with 15 and 30 K grey- 
body components. We discuss this further in Sect. 7, in 
terms of the geometry and opacity of the disk. We also 
note that in passing from the centre to the outer regions, 
progressively bigger grains exhibit stochastic heating, as 
expected due to the lower energy density of the radiation 
field. For the same location in the disk, graphite grains 
undergo larger fluctuations than the silicate grains. 

Given the temperature probability function P(a,T), 
the IR flux is given by 

Fv = j2 J N(a)dair a 2 Q &hs (v,a) 

x / B,{T)P(a,T)dT (16) 
Jo 

where B v is the Planck function, N(a) is the grain size 
distribution given by Eq. (1), and d is the distance to the 
galaxy. 

The IR emission was calculated in each point of the 
galaxy and then the total emission was obtained by in- 
tegrating the emission over the galaxy. In Fig. 3 we give 



two examples of IR spectra for two extreme locations in 
the plane of the galactic disk, namely in the center of the 
galaxy (R = kpc, z = kpc) and at the edge of the disk 
(R = 15 kpc, z = Okpc). At the edge of the disk the in- 
frared emission is colder, as expected due to the decrease 
in the dust temperature, though the increased importance 
of stochastic heating in the outer regions tends to flatten 
the SED on the Wien side. 

3. Results for the "standard model" 

The FIR emission for our "standard model" was calcu- 
lated for a grid of SFR, and the different solutions were 
compared with the observations. To extract the FIR lu- 
minosity of the galactic dust we integrated the SED of 
NGC891 between 40 and 1000 /im, with the data taken 
from Alton et al. (1998). The SED was parameterised by 
Alton et al. (1998) by a fitting two-temperature greybody 
curve (with a wavelength dependence of grain emissivity 
in the FIR (5 = 2) to the observed fluxes from IRAS and 
SCUBA at 60, 100, 450 and 850 /im. This yielded best-fit 
"temperatures" of 30 and 15 K to the observed spectrum, 
and an integrated luminosity between 40 and 1000 of 
5.81 x 10 36 W. This output luminosity can be accommo- 
dated in our "standard model" with a solution for a SFR 
= 7.5M Q /yr. The predicted FIR SED from the model is 
given in Fig. 4a, together with the UV-optical-NIR intrin- 
sic emitted stellar radiation (as would have been observed 
in the absence of dust). The IRAS and SCUBA data (dia- 
monds) from Alton et al. (1998) are also given for compar- 
ison. The inferred SFR is an upper limit, since we assumed 
that all the non-ionising UV escapes into the diffuse com- 
ponent. Indeed, in the absence of local absorption, all the 
lines of sight from each source would contribute primary 
photons to the diffuse interstellar radiation field. In prac- 
tice, some fraction of the lines of sight - call this fraction 
F - from each source will be opaque to local absorption, 
giving rise to discrete IR sources embedded in a diffuse IR 
emitting disk. Our solution was thus calculated under the 
extreme assumption that F = for the non-ionising UV 
and F = 1 for the ionising UV. While requiring a solution 
with a large SFR, the predicted SED is slightly warmer 
than the observed SED at 100 /im. However, the main dif- 
ficulty of the model is that it underestimates by a factor 
of 2 to 4 the observed sub-mm fluxes. 

4. Results for the "standard model" 
supplemented by localised sources 

A more realistic solution is to consider that some frac- 
tion F of the non-ionising UV is locally absorbed in star- 
forming complexes (e.g., HII regions). Such a solution re- 
quires two free parameters (F and SFR) and can be ac- 
commodated with a lower SFR, due to the higher prob- 
ability of absorption of the non-ionising UV photons. To 
include the contribution of HII regions we used a template 




Fig. 4. a) The SED from our "standard model", for a current SFR = 7.5M /yr. The UV-optical-NIR SED plotted in 
the left hand panel represents the intrinsic emitted stellar radiation (as would have been observed in the absence of dust) 
while the FIR-sub-mm curve plotted in the right hand panel is the predicted SED for the re-radiated dust emission. 
The SFR has been adjusted such that the predicted FIR model luminosity equals the observed one. The observed 
FIR fluxes from Alton et al. (1998) are given as diamonds, b) The SED from our "standard model" supplemented 
by localised sources, for a SFR = 3.5M /yr and F = 0.28. The UV-optical-NIR SED plotted in the left hand panel 
represents the intrinsic emitted stellar radiation while the FIR-sub-mm curve (solid line) plotted in the right hand 
panel is the predicted SED for the re-radiated dust emission. In the right hand panel we also show the contribution of 
the diffuse component (dashed-line) and of the HII component (dotted line) to the total predicted FIR emission. The 
observed FIR fluxes from Alton et al. (1998) are again given as diamonds. 



SED of such forming complexes. We utilised the galac- 
tic UC HII region G45.12+0.13, with FIR observations at 
60 and 100 /jm from the IRAS Point Source Catalog and 
1300 fim measurements from Chini et al. (1986). We fitted 
the observed data with a grcybody curve (with a wave- 
length dependence of grain emissivity in the FIR [3 = 2) 
and used it as a template SED for the compact source 
component for our modelling. We do not attempt here 
to include in this template potential cold emission com- 
ponents that might be expected from "parent" molecular 
clouds in juxstaposition to their "offspring" HII regions. 

The best solution was obtained for a SFR =3.5 M Q /yr 
and F = 0.28. The predicted FIR SED is given in Fig. 4b, 
again together with the UV-optical-NIR intrinsic stellar 
radiation. In Fig. 4b we also give the contribution of the 
diffuse component (dashed line) and of the HII component 
(dotted line) to the the total predicted FIR SED (solid 
line). The integrated emitted luminosity predicted by the 
model is L = 4.79 x 10 36 W, which accounts for 82% of 
the total observed emitted FIR luminosity. Of this, the 
diffuse component contributes 46% of the total observed 
dust luminosity and the localised sources 36%. However, 
the main difficulty of the model is again that the solution 
fails to reproduce the observations longwards of 100 ^m. 



5. The missing FIR/sub-mm component of our 
"standard model" 

We have shown that our "standard model" fails to repro- 
duce the observed SED of NGC891 in the sub- mm spec- 
tral range. The more realistic case, where the contribution 
from HII regions is added to the diffuse component, is also 
not able to predict the observed FIR/sub-mm SED. How- 
ever, this may be, to some extent, because we did not 
include a cold dust emission component that might be ex- 
pected to arise from associated parent molecular clouds. 

Below we consider four more ways that might account 
for this discrepancy. Firstly, we discuss the effect of alter- 
ing the geometry of the large-scale distribution of stars 
and dust. Alternatively, we consider extraplanar emis- 
sion. Thirdly, the predicted apparent sub- mm/optical ra- 
tio could be increased by including more grains in the 
disk, but distributed in clumps, so as not to effect the 
optical extinction properties of the diffuse disk. Lastly, a 
more radical solution is presented, namely a model with 
two dust-disk components, which is a moderately optically 
thick solution. 

5.1. The effect of altering the geometry of the large-scale 
distribution of stars and dust 



The diffuse IR emission predicted by our model was de- 
rived for an exponential disk of stars and dust, adopted 
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from Xilouris et al. (1999). Galactic disks are however 
known to be quite complex systems, where the large-scale 
distribution of stars and dust presents inhomogeneities in 
the form of spiral structures. Adding logarithmic spiral 
arms as a perturbation on the exponential disks for stars 
and dust would alter the solutions of the radiative transfer, 
and thus of the calculated diffuse IR emission. However, 
Misiriotis et al. (2000b) showed that plain exponential disk 
models give a very accurate description for galactic disks 
seen edge-on, with only small deviations in parameter val- 
ues from the real ones (typically a few percent). Thus, 
we conclude that the geometry adopted here cannot be 
responsible to account for the discrepancy between our 
model predictions and the observations. 

5.2. Extraplanar dust? 

Observations of external galaxies have revealed extensive 
thickened layers of ionised gas traced by its Hen emission 
in several edge-on spiral galaxies (e.g., Rand et al. 1990, 
1992; Dettmar 1990; Pildis et al. 1994; Rand 1996, Fer- 
guson et al. 1996; Hoopes et al. 1999). Amongst the most 
spectacular examples of extraplanar diffuse ionising gas 
(DIG) is that seen in NGC891. The recent observations of 
Hoopes et al. (1999) showed that the DIG layer in NGC891 
extends out to at least 5 kpc from the plane, and possibly 
as far as 7 kpc, and has the brightest and largest DIG 
layer known. 

In addition to detecting the high z DIG studied in 
earlier work, Howk & Savage (2000) also detected individ- 
ual dust-bearing clouds observable to heights z ~ 2 kpc 
from the midplane. Moreover, they detected the presence 
of HII regions at large distances from the midplane (0.6- 
2 kpc), which suggests that on-going star formation may 
be present in some of the dense, high-z clouds. 

In principle both the dust-bearing clouds (z <~ 2 kpc) 
seen in absorption as well as the dust that may be as- 
sociated with the DIG (z ~ 5 kpc) could contribute to 
the FIR/sub-mm emission of the galaxy. The possibility 
of the extraplanar dust clouds emitting in the sub-mm 
was recently investigated by Alton et al. (2000b). From 
the upper limits in the total amount of extraplanar dust 
these authors found that less than 5% of galactic dust ex- 
ists outside the galactic disk, if the dust grains are not 
colder than 17 K, or 9%, if the dust temperature is 10 K. 
This small percentage cannot account for our discrepancy 
between the model and the observations. 

We also consider it unlikely that any diffuse dust com- 
ponent associated with the DIG (z ~ 5 kpc) could pro- 
vide a substantial fraction of the sub-mm emission. In 
that case, if the extraplanar emission were powered by 
the absorption of disk photons, the optical depth of the 
DIG would have to be comparable to that of the disk, 
which would be a very extreme scenario, probably also in 
conflict with the constraints on submillimetcr extraplanar 
emission. 



5.3. The case of hidden dust in clumps 

One possible explanation to account for the missing dust 
component could be that Xilouris et al. (1999) underes- 
timated the dust content of the galaxy because of the 
adopted diffuse distribution for the dust. Since the galaxy 
is optically thin, any additional dust component would 
have to be distributed in very dense clumps instead of be- 
ing diffuse. Only then would a significant fraction of the 
lines of sight avoid the additional dust, so as not to affect 
the derived optical-NIR optical depth. 

Such dust clumps could either be a component of star- 
formation regions, or could have no associated sources - we 
refer to the latter hypothesis as the "quiescent clumps". 
The quiescent clumps must be optically thick to the diffuse 
UV/optical radiation field in the disk to have an impact 
on the predicted submillimcter emission. They would radi- 
ate at predominantly longer wavelengths than the diffuse 
disk emission, in the FIR/sub-mm spectral range. One 
can speculate that such optically thick "quiescent clouds" 
could be physically identified with partially or wholly col- 
lapsed clouds that, for lack of a trigger, have not (yet) 
begun to form stars. However, due to the lack of intrin- 
sic sources, a very substantial mass of dust would have to 
be associated with the quiescent clouds to account for the 
submillimetcr emission missing from our model prediction. 

In reality there must be also dark clouds associated 
with star-forming complexes. In the Milky Way HII re- 
gions around newly born massive stars are commonly seen 
in juxtaposition to parent molecular clouds (e.g., M17). 
This is thought to be a consequence of the fragmenta- 
tion of the clouds due to mechanical energy input from 
the winds of the massive stars. Thus, warm dust emis- 
sion from cloud surfaces directly illuminated by massive 
stars can be seen along a fraction of the lines of sight, to- 
gether with cold sub-mm dust emission from the interior 
of the associated optically thick cloud fragments. These 
dense clumps (with or without associated sources) may 
contribute to the sub-mm emission, and thus supplement 
the contribution of compact HII regions. 

However, the question remains whether it is possible 
to place such a substantial amount of additional dust in 
clumps without affecting the opacity of the galaxy. Kuhin- 
ski et al. (1998) modeled a sample of highly inclined 
galaxies using both smooth and clumpy dust and con- 
cluded that dumpiness does not effect dramatically the 
opacity of highly inclined galaxies. On the other hand, 
Bianchi et al. (1999) adopt different parameterisation for 
clumps and they report that dumpiness does significantly 
effect the opacity of edge-on galaxies. It is beyond the 
scope of this paper to calculate a quantitative solution for 
a clumpy distribution, but we qualitatively consider the 
clumpy scenario, as quiescent clumps or associated with 
star-formation regions - to be a possibility to account for 
the missing FIR-sub-mm component. 
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5.4- A two-dust-disk model 

We have seen that the observed sub-mm SED requires 
more dust in the galaxy than predicted by our "standard 
model". We have already discussed the possibility of in- 
cluding this extra dust in chimps. An alternative solution 
is obtained if it is postulated that this additional dust is 
confined in a diffuse thin disk, in which the young stel- 
lar population is embedded. (A thin young stellar disk, 
but with no associated dust, was considered in Xilouris et 
al. 1998.) If this second dust-disk component has a scale- 
height comparable with the 90 pc scaleheight of molecular 
gas clouds in the Milky Way, it would not have been seen 
in the optical/NIR analysis of Xilouris et al. (1999), being 
totally obscured by the optically visible disk of dust. In- 
deed, since the galaxy is almost edge-on, adding a second 
thin disk of dust will probably not change the parameters 
for the intrinsic distributions of stars and dust derived by 
Xilouris ct al. (1999). However, adding a second disk of 
dust will change the optical depth of the disk, transform- 
ing our solution into a moderately optically thick solution 
(4 = 3.1). 



100.0 r 
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Fig. 5. The predicted FIR emission from a two-dust-disk 
model, for a SFR = 3.8M /yr and F = 0.22. The total 
predicted SED is given with solid line, the diffuse compo- 
nent with dashed line and the HII component with dotted 
line. The observed fluxes from Alton et al. (1998) are given 
as diamonds. 

We have calculated the energy density of the radia- 
tion field from UV to NIR for such a two-dust-disk model, 
and subsequently the infrared emission. The best solu- 
tion was obtained for a SFR = 3.8M /yr, F = 0.22 and 
Mdust — 7 x 10 7 M Q in the second disk of dust. The cor- 
responding non-ionising UV luminosity for this SFR is 
~ 8.2 x 10 36 W. This could be provided, for example, by 
a population of B5 stars (T e g = 15500 K) with a space 
density in the centre of the galaxy of 7.5 x 10~ 4 pc~ 3 . 
The predicted combined SED from the two diffuse dust- 
disk components is shown by the dashed line in Fig. 5. 
By adding the further contribution from localised sources 
(dotted line), as discussed in Sect. 4, we were able to fit 



the SED of NGC891, from the FIR to the sub-mm. The lu- 
minosity of the diffuse component is 4.07 x 10 36 W, which 
accounts for 69% of the observed FIR luminosity, and the 
luminosity of the HII component is 1.82 x 10 36 W, mak- 
ing up the remaining 31% of the FIR luminosity This 
model can thus successfully fit the shape of the SED. In 
the next section we will show that the predictions of the 
two-dust-disk model are in agreement with the observed 
radial profile at 850 /mi. Nevertheless, further modelling 
of face-on galaxies will be required to test the validity 
of the two-dust-disk model. Ultimately, the two-dust-disk 
hypothesis will be directly tested observationally using the 
new generation of sub-mm interferometers, which will re- 
solve edge-on disks. 

6. The radial profiles 
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Fig. 6. The averaged radial profile at 850 /mi from the 
two-dust-disk model, plotted with the solid line. The pro- 
file is averaged over a bin width of 36", for a sampling of 
3" and for a beam width of 16", in the same way as the 
observed averaged radial profile from Alton et al. (2000a) 
(plotted with diamonds). 

Since we have calculated the 3 dimensional FIR radia- 
tion field in NGC891, for both our "standard model" and 
the two-dust-disk model, it is straightforward to integrate 
along the line of sight and thus produce a 2-dimensional 
map of the galaxy at different wavelengths, and then ra- 
dial profiles. Here we consider only the diffuse component 
of the IR emission. Because observed radial profiles were 
derived by Alton et al. (2000a) using the SCUBA obser- 
vations at 850 /mi, we first attempt to calculate the radial 
profiles at this wavelength and compare it with the ob- 
servations. The calculations were done at high resolution 
and afterwards smoothed to 16" resolution. Since the ob- 
served radial profile at 850 /mi from Alton et al. (2000a) 
were derived as averages over a bin width of 36", and with 
a sampling of 3", we have also calculated an averaged ra- 
dial profile in the same way as Alton et al. (2000a) did, 
and directly compare it with the observations. We have 
found that in the case of the two dust disk model there is 
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Fig. 7. a) The predicted diffuse FIR radial profile from our "standard model" at 850 /mi, for SFR = 3.5M Q /yr. The 
radial distance is given in arcsec, with 100" ~ 5 kpc). b) The predicted colour profiles at different wavelengths, for the 
same model like in a). The ratios were multiplied by the given factors to provide the same colours at an angular radius 
of 50". The colour profiles show that for shorter wavelengths the effect of decreasing the dust temperature towards 
larger radii become more and more pronounced. The colour profile F60/F850 has the steepest gradient, with a very 
hot component within the inner 50 arcsec radius, followed by the increase of the cold dust at the larger galactocentric 
distances, c) The predicted diffuse FIR radial profile from our two-dust-disk model, at 850 /im, for SFR = 3.8M©/yr. 
d) The colour profiles for the two-dust-disk model. Again, the ratios were multiplied by the given factors to provide 
the same colours at an angular radius of 20". 



a very good agreement between the model predictions and 
the observations and this comparison is illustrated in Fig. 
6. The predicted radial profile can be traced out to 300 
arcsec radius (15 kpc), as also detected by the SCUBA. 

In Figs. 7a and 7c we give our prediction for the ra- 
dial profile of NGC891 at 850 /im, this time as it would 
be observed within a beam of 16", and without any aver- 
aging. This is shown for our "standard model" and the 
two-dust-disk model, respectively. In Fig. 7b. d we give 
the radial colour profiles at the wavelengths used by the 
IRAS, ISO and SCUBA, normalised to the radial profile 
at 850 /xm. The strong colour gradients, seen in partic- 
ular between the 60 or 100 (im and the sub-mm regime, 
are a consequence of the strong gradients in the radia- 



tion field calculated in our finite disk model. It is inter- 
esting to note that the stronger contribution of stochastic 
emission from small grains (relative to big grains) in the 
outer disk at 60 /im, does not strongly influence the gra- 
dients. Potentially, therefore, large-scale colour gradients 
in galaxy disks, could act as a diagnostic for the relative 
contribution of the diffuse emission and compact sources, 
as the latter would not be expected to have IR spectra 
strongly dependent on position in the disk. As expected, 
the 450/850 colour shows only a shallow gradient, as both 
these wavelengths are displaced well longwards of the spec- 
tral peak, and thus depend only approximately linearly on 
grain temperature. 



12 



C.C.Popescu et al. 



7. Discussion: the contribution of different stellar 
populations in heating the dust in NGC891 

The relative contribution of optical and UV photons in 
heating the dust has been a long standing question in the 
literature. Since we have a detailed calculation of the ab- 
sorbed energy over the whole spectral range and at each 
position in the galaxy, we can directly calculate which part 
of the emitted FIR luminosity from each volume element 
of the galaxy is due to the optical and NIR photons, and 
which part is due to the UV photons. The IR emission 
from any given volume clement arises from a common lo- 
cal grain population, so that the IR colours from that vol- 
ume element will be the same for the emitted IR photons 
regardless of whether the absorbed energy had an optical 
or UV origin. However the relative fraction of energy ab- 
sorbed from the UV and optical light varies with position 
in the galaxy, as does the dust temperature, so that the 
volume-integrated IR spectrum representing re-radiated 
optical light will differ from that of the re-radiated UV 
light. In this way, some information on the stellar popu- 
lations heating the grains is preserved in the IR domain, 
in the form of the observed IR colours. Volume-integrated 
IR spectral components arising from re-radiated optical 
and UV light are presented in Fig. 8a, b for the "standard 
model" and in Fig. 8c, d for the two-dust-disk model, re- 
spectively. The overall contribution of the diffuse optical 
and NIR radiation (4000-22000 A) to the total FIR lumi- 
nosity is 1.07 x 10 36 W, or 22% for the "standard model". 
Since the predicted luminosity of this model is lower than 
the observed one, it is also meaningful to give the contri- 
bution of the diffuse optical radiation to the observed FIR 
luminosity, which is 18%. This component of the FIR SED 
is plotted with dashed line in Fig. 8a, b. The correspond- 
ing optically generated IR luminosity for the two-dust- 
disk model is 1.85 x 10 36 W, or 31%, higher than in the 
"standard model" , largely because the galaxy is becoming 
optically thick in the optical range. This component of the 
FIR SED is again plotted with dashcd-linc in Fig. 8c, d. 
The contribution of different stellar components to the to- 
tal FIR luminosity for the two dust models is summarised 
in Table 2. Table 3 gives the contribution of different stel- 
lar components to the FIR wavelength commonly used by 
IRAS, ISO and SCUBA. 

We note that, regardless of the optical thickness of the 
adopted model, the diffuse optical radiation field makes 
only a relatively small contribution to the total emitted 
dust luminosity. This is in qualitative agreement with var- 
ious statistical inferences linking FIR emission with young 
stellar populations, in particular the FIR-radio correlation 
which, also after normalisation to indicators of galactic 
mass, such as K-band flux, remains rather tight (Xu et al. 
1994). Perhaps more surprising, at first glance, is the pre- 
dicted predominance of UV-powered grain emission in the 
sub-mm range in both our dust models (Fig. 8b, d.). How- 
ever, our analysis has demonstrated that the FIR colours 



Table 2. The contribution of the diffuse optical/NIR ra- 
diation L°Pj R (4000 -22000 A), diffuse UV radiation LpJ R 
(912-4000 A) and UV photons locally absorbed in HII re- 
gions Lpj R to the total predicted FIR luminosity L^ IR l ■ 
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Table 3. The relative contribution of the diffuse opti- 
cal/NIR radiation f°f R (4000 - 22000 A), diffuse UV ra- 
diation f%Y R (912 - 4000 A) and UV photons locally ab- 
sorbed in HII regions fp/ R to various FIR wavelengths. 
The wavelengths considered here were those commonly 
used by IRAS, ISO and SCUBA, namely 60, 100, 170, 
450, and 850 fim. 
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have to be fundamentally interpreted also in geometrical 
terms, rather than simply as separate temperature com- 
ponents. Thus, the predicted increasing predominance of 
re-radiated diffuse UV photons over re-radiated diffuse op- 
tical photons going from the FIR into the sub-mm regime 
is due to an increasing fraction of the emission arising from 
more optically thin regions, where the ratio of optical to 
UV absorption is largely controlled by the relative optical 
and UV grain cmissivities. This particularly applies in the 
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Fig. 8. a) The contribution of different stellar components to the FIR emission, as a function of wavelength, for our 
"standard model". The legend is as follows: dashed- line: contribution of diffuse optical radiation (4000 — 22000 A), 
dashed-dotted line: contribution of diffuse UV radiation (912 — 4000 A), dotted-line: contribution of HII component, 
solid line: total predicted FIR SED. The observed fluxes from Alton et al. (1998) are given as diamonds, b) The 
same as in case a), but here the different contributions are given as the fraction f of the FIR emission produced by a 
certain component to the total FIR predicted emission, c) and d) are the same as a) and b), respectively, but for the 
two-dust-disk model. 



the outer disk, which, in addition to being (relative to the 
inner disk) optically thin, has systematically lower grain 
heating and temperature. 



8. Summary and Conclusions 

In this paper we have modeled the edge-on spiral galaxy 
NGC891 from the UV to the FIR and sub-mm with 
the purpose of understanding what is the origin of the 
FIR/sub-mm emission. The strategy of the paper was to 



use a model for the intrinsic distribution of older stars and 
associated dust as derived from fitting the optical and NIR 
images of NGC891 by Xilouris et al. (1999), which, sup- 
plemented by a distribution of newly- form stars, should 
predict the diffuse FIR emission. By comparing the pre- 
diction of the "standard model" with the observations we 
tried to identify whether this model can reproduce the 
observed dust emission. 

To calculate the FIR emission in the galaxy we needed 
to have a dust model and also the energy density of the 
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radiation field at each position in the galaxy. In describing 
the properties of the grains we used the extinction efficien- 
cies from Laor & Draine (1993) for spherical "astronomi- 
cal grains", for a mixture of graphites and silicates of the 
MRN grain size distribution. For small grains not in equi- 
librium with the radiation it was necessary to model the 
stochastic emission. For this calculation we adopted the 
heat capacities from Guhathakurta & Draine (1989) for 
silicate grains and from Dwek (1986) for graphite grains. 
Our theoretical dust model was found to be consistent 
with the extinction coefficients in the optical and NIR 
spectral range, derived independently by Xilouris et al. 
(1999) for NGC891. The similarity between the extinction 
coefficients derived from our theoretical model and those 
derived as fitting parameters assure us that we can use 
the model for the distribution of old stars and associated 
dust of Xilouris et al. (1999). 

The calculation of the energy density in the galaxy was 
done using the solution of the radiation transfer equation, 
based on the method of Kylafis & Bahcall (1987). This cal- 
culation was done separately for the optical/NIR and UV 
radiation field. For the distribution of old and intermediate 
age stars we used a smooth exponential disk and a de Vau- 
couleurs bulge, with the parameters derived from Xilouris 
et al. (1999). For the distribution of newly-formed stars 
we used a smooth exponential disk, with fixed parameters, 
e.g., with the same scalelengths as that of the stars radi- 
ating in the B band, and with a scaleheight of 90 pc. The 
amplitude of the radiation emitted by the young stars (of 
the ionising and non-ionising UV radiation) was consid- 
ered as a free parameter. Using the population synthesis 
models of Mateu & Bruzual (2000) we parameterised the 
UV radiation in terms of the recent SFR. Another free 
parameter was the fraction F between the UV photons 
which are locally absorbed within star-forming complexes 
and those which participate in the diffuse ISRF. 

The predicted SED from our "standard model" was 
compared with the observed one. In the extreme case that 
F = 0, the SFR can be taken so as to exactly reproduce 
the total bolomctric output in infrared. This is a solution 
with SFR = 7.5M /yr. While requiring a high SFR, this 
solution fails to reproduce the sub-mm part of the spec- 
trum. A solution with F > 0, which includes also the con- 
tribution of star- forming complexes, underestimates even 
more the emission in the sub-mm spectral range. The fail- 
ure of our "standard model" to fit the observed SED of 
NGC891 indicates, most probably, that there is missing 
dust, which could not have been detected by Xilouris et al. 
(1999). Two scenarios are proposed to explain this missing 
dust component, namely a scenario in which the missing 
dust is in the form of very small optically thick clumps, 
such that they do not affect the optical extinction derived 
by Xilouris et al. (1999), or in the form of a second diffuse 
thin disk, where the young stellar population is embedded. 
For the second case we give a quantitative description of 
the model, and we make calculations for the FIR output. 



We have found that the two-dust-disk model reproduces 
very well the observed SED, as well as the radial profile 
at 850 ^m. This solution requires a SFR = 3.5M Q /yr and 
F = 0.28. This value of the SFR fits better with the pre- 
conception that NGC891 is a relatively quiescent normal 
galaxy. 

Potentially, a further possibility exists to explain the 
shortfall in sub-mm emission predicted by our "standard 
model" . A higher level of sub-mm emission could be at- 
tained from a given grain mass if the sub-mm grain emis- 
sivity were higher than assumed in our dust model. This 
approach was followed by Alton et al. (2000a). However, 
for a given grain heating the sub-mm/FIR colour of the 
grain emission will become colder, in order to preserve 
the energy balance between the absorbed UV/optical and 
the radiated emission. Thus, the peak of the SED would 
be shifted to longer wavelengths, so in this scenario there 
would be a need to increase the SFR to account for the 
observed fluxes in the 60-100 /im range. 

We emphasise that, on the basis of currently avail- 
able observational evidence, it is difficult to distinguish 
between the two dust-disk-model and the clumpy scenario. 
As present day sub-mm and FIR telescopes generally have 
insufficient angular resolution to resolve the postulated 
geometrical emission components, the best way of distin- 
guishing between them might be to compare the statistics 
of quantities derived from the UV/optical/FIR/sub-mm 
SEDs. Such statistics could be, for example, the derived 
star-formation rates from a sample of edge-on and face-on 
systems with comparable morphological types. With the 
advent of the future generations of FIR telescopes with 
arcsec resolution we will be able to directly distinguish 
between diffuse disk and cloud scenarios. Then such tech- 
niques will be needed to understand young galaxies at 
cosmological distances. 
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